Magnetic field gradients in solar wind plasma and geophysics periods 
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Using recent data obtained by Advanced Composition Explorer (ACE) the pumping scale of 
the magnetic field gradients of the solar wind plasma has been calculated. This pumping scale is 
found to be equal to 24h ± 2h. The ACE spacecraft orbits at the LI libration point which is a 
point of Earth-Sun gravitational equilibrium about 1.5 million km from Earth. Since the Earth's 
magnetosphere extends into the vacuum of space from approximately 80 to 60,000 kilometers on the 
side toward the Sun the pumping scale cannot be a consequence of the 24h-period of the Earth's 
rotation. Vise versa, a speculation is suggested that for the very long time of the coexistence of 
Earth and of the solar wind the weak interaction between the solar wind and Earth could lead to 
stochastic synchronization between the Earth's rotation and the pumping scale of the solar wind 
magnetic field gradients. This synchronization could transform an original period of the Earth's 
rotation to the period close to the pumping scale of the solar wind magnetic field gradients. 
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INTRODUCTION 



Magnetic field in solar wind plasma is actively studied 
in the last years both theoretically and experimentally 
(see, for instance, Properties of this field 

gradients are of especial interest because of strong non- 
homogeneity of the field. Inferring universal properties of 
the magnetic field is a difficult task due to superposition 
of the strong non-homogeneity and global anisotropy. 
Even in the inertial range of scales the non-homogeneity 
and anisotropy affect behavior of different components of 
the magnetic field and its gradients. It can be shown 
however, that magnitude of the magnetic field B = \J Bf 
does exhibit certain universal properties in inertial range 
of scales. Moreover, we will show in present paper that 
magnitude of gradients of the solar wind magnetic field 
exhibits certain universal properties as well. These 
universal properties have substantial geophysical conse- 
quences, which we discuss in the last section of the paper. 

In magnetohydrodynamics (MHD) the magnetic field 
fluctuation B dynamics is described by equation 

— = V x (v x B) + r;V 2 B. (1) 

Here, v is the turbulence velocity and the r\ is magnetic 
diffusivity. Equation (1) can be regarded as a vector 
analogue of the advection-diffusion equation 

^-(v-v)« + j)v 2 e (2) 

at 

for the evolution of a passive scalar subject to molecular 
diffusivity D. Aside from the fact that B is a vector and 
a scalar, the equations are different also because v in 
Eq. (1) can be affected quite readily by the feedback of 
the magnetic field B. Our interest here is to explore the 
extent of similarities, despite these obvious differences, in 
the inertial (Batchelor) range statistics of the magnetic 




FIG. 1: Energy spectrum of the magnitude B = y/Bf of the 
magnetic field B in the solar wind plasma, as measured by 
the ACE magnetometers in the nanoTesla range for the year 
1998 (4 min average). 



and the passive scalar fields (see, for instance, Refs. [151 ]- 
[l8| and the papers cited there). 

Solar wind is an excellent natural "laboratory" for the 
MHD problem. It is known that the statistical prop- 
erties of velocity fluctuations in the solar wind are re- 
markably similar to those observed in fluid turbulence 
[l9| . It is also known that the plasma power spectra of 
the magnetic field and velocity fluctuations often contain 
an "inertial" range with a slope of approximately —5/3 
(for reviews, see [l!|,[2(|). The approximately —5/3 
power-law is especially common for magnitude fluctua- 
tions B = \/Bf (the summation over repeated indexes 
is assumed) of the magnetic field, as one can see in Fig. 
1. For computing this spectrum, we have used the data 
obtained from Advanced Composition Explorer (ACE) 
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satellite magnetometers for the year 1998. In this pe- 
riod, the sun was quiet and the data are statistically sta- 
ble. The range of scales for which the "—5/3" power 
holds is taken to be the inertial range; the smaller scales 
are obliterated because of the instrument resolution and 
the truncation at the large-scale end is governed by the 
record length chosen for Fourier transforming. 

The nature of the spectrum for each individual com- 
ponent of the magnetic field is more variable from one 
component of B to another, and from one situation to an- 
other, perhaps because of large anisotropics in the mag- 
netic field B, but the result for the magnitude of B seems 
more robust. Scaling spectrum with the —5/3 slope (Fig. 
1) is quite typical of that observed for passive scalar fluc- 
tuations in fully developed three-dimensional fluid tur- 
bulence (the so-called Corrsin-Obukhov spectrum [22I]). 
Spurred by this similarity, we were motivated to explore 
further the properties of the magnitude B and its gradi- 
ents, and compare them with those of the passive scalar. 



STATISTICAL PROPERTIES OF THE 
MAGNITUDE OF THE MAGNETIC FIELD 

More detailed statistical information is provided by the 
structure functions scaling 



where 



AB T = B(t + r) - B(t). 



(3) 



(4) 



The exponent £2 is directly related to the spectral ex- 
ponent (in our case £2 ss 5/3 — 1 = 2/3 [22j). If the 
dependence of Cp on p is nonlinear, it is well-known that 
one has to deal with intermittency. 

Figure 2 shows the scaling of structure functions for 
the solar wind data. Slopes of the straight-line fits in 
the apparently scaling region provide us the scaling ex- 
ponents these are shown in Fig. 3 as circles. Trian- 
gles in the figure indicate experimental values obtained 
for temperature fluctuations in the atmosphere [23j . The 
other experimental data [25j],[26| are in agreement with 
each other to better than 5%. The * symbols are for the 
passive scalar field obtained by numerically solving the 
advection-diffusion in three-dimensional turbulence [24 1 . 



It is clear that the exponents for the passive scalar data 
are in essential agreement with those for the magnitude 
fluctuations of the magnetic field. 

One can, in fact, analyze the solar wind data somewhat 
differently using the notion of the extended self-similarity 
(ESS). Since, empirically, the fourth order exponent is 
quite closely equal to 1 for the magnetic field magnitude, 
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FIG. 2: Structure functions of magnetic field magnitude in the 
solar wind plasma as measured by the ACE magnetometers 
in nanoTesla for the year 1998 (4 min averages). The straight 
lines (the best fits) are drawn to indicate the scaling law (3) 
in the inertial range. 
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FIG. 3: Scaling exponents (3) calculated for B in the solar 
wind (circles) and for the passive scalar in the atmospheric 
turbulence (triangles, [23j), and in the direct numerical sim- 
ulation of 3D fluid turbulence (stars, [241 ]) 



we can extend the scaling range (and consequently im- 
prove the confidence with which those exponents are de- 
termined) by redefining them as 



(|AS T H ~ (|AS r 



4\C P 



(6) 



(|AB T | 4 )~r, 



(5) 



Figure 4 shows the ESS dependence (6). The slopes of 
the best-fit straight lines in this figure provide us with 
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FIG. 4: Extended self-similarity (ESS) of the magnetic field 
magnitude in the solar wind plasma. The straight lines (the 
best fit) are drawn to indicate the ESS (6). 



the magnetic field are essentially the same. This is both 
surprising and thought-provoking, and needs to be un- 
derstood further. To this end, let us return to Eq. (1) 



and specialize 28], for simplicity, to the case of incom- 
pressible moving medium (V • v = 0). Equation (1) can 
then be rewritten as 



<9B 

~dt 



-(v V)B + (B- V)v- 



ryV 2 B. 



(7) 



Let us now consider the equation for the magnitude B of 
the magnetic fluctuations given by B = Bn, where n is 
the unit vector with its direction along B: m = Bi/B. 
Multiplying both sides of Eq. (7) by the vector n and 
taking into account that nf = 1 we obtain 



dB 
~dt 



-(v- V)B + r/V 2 B + \B, 



(8) 



in which the "friction-stretching" (or the production) co- 
efficient A in the last term has the form 



A 



dvj 
' dxj 



drii 



(9) 
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FIG. 5: The same as in Fig. 3 but using the ESS method (6) 
for B. 



the ESS scaling exponents £ p , which are shown in Fig. 5 
as circles. The other symbols have remained unchanged 
from Fig. 3. The shift of the exponents £ p in comparison 
to those from ordinary self-similarity is about 4%, but 
the scaling interval for ESS is considerably larger. This 
increased scaling range is well-known in other contexts 

The results shown in Figs. 3 and 5 suggest that at 
least up to the level of the fourth-order the scaling ex- 
ponents for the passive scalars and for the magnitude of 



with the indexes i and j representing the space coordi- 
nates, and the summation over repeated indexes is as- 
sumed. The first term on the right hand side of Eq. (9) 
is crucial for any dynamo effect. 

If the statistical behaviors of 9 and B are to be similar, 
as suggested by Figs. 1, 3 and 5, we should be able to 
observe the underlying similarity between Eqs. (2) and 
(8). There is a major difference corresponding the pres- 
ence in Eq. (8) of the production term XB. However, 
given the empirical indications that B and 9 arc similar 
in the inertial range, it is appropriate to look for circum- 
stances under which the A term in Eq. (8) may be small. 
The second term in A is assured to be small because the 
smallness of the magnetic diffusivity 77, but difficulties 
may arise from the first term on the right hand side of 
Eq. (9). 

To eliminate directional dependencies in Eq. (8), let us 
make the following conditional average of that equation. 
That is, fix the magnitude B in the vector field B = Bn 
while performing the average over all realizations of the 
direction vector field n permitted by the vector equation 
(7). Let us denote this ensemble average as (...) n . From 
the definition, this averaging procedure does not affect 
B itself, but affects the velocity field v and the "friction- 
stretching" coefficient A in Eq. (8). We thus obtain 



dB 



V)B 



f]V 2 B + (\) n B. (10) 



It is worth emphasizing that the solutions of the original 
equation (7) satisfy Eqs. (8) and (10), but not all possi- 
ble formal solutions of the Eqs. (8) and (10) satisfy Eq. 
(7); similarly, not all formal solutions of Eq. (10) sat- 
isfy Eq. (8) while all solutions of Eq. (8) do satisfy Eq. 



4 



(10). Restricting comments to the relationship between 
Eqs. (8) and (10), the solutions of the two equations are 
the same only if the initial conditions are the same and 
if realizations of (v) n and of (A) n , related to these ini- 
tial conditions by the conditional average procedure, are 
obtained from solutions applicable to Eq. (8). 

Returning now to Eq. (10), the conditionally averaged 
velocity field (v) n may posses statistical properties that 
are different from those of the original velocity field v, 
and there can be circumstances under which (A) n = 0, 
or small. If so, the similarity between Eqs. (2) and (10) 
(and, consequently, Eq. (8)) can be the basis for the sim- 
ilarity in statistical properties of their solutions. There- 
fore, finding conditions under which (A) n = 0, or small, 
seems to be a useful exercise. 

It is, however, difficult to guess a priori when (A) n is 
negligible, because there is no small parameter for the 
stretching part of A. Therefore, let us consider a generic 
set of conditions, presumably for the inertial range, which 
can result in {niUjdvi/ dxj) n = 0. This can be a combi- 
nation of isotropy, which yields 

(n i n j ) I1 = (i^j) 

and 

(n?>„ = (n|> n = <n§) n , (11) 
and statistical independence 

{riiTlj^n = (niTlj} a (tp} a , (12) 

where tp = dvk/dxi for arbitrary k and I. 

We should emphasize that the conditional average in- 
dicated by (. . .)„ and the global average indicated by 
(. . .) are quite different; because of this, the quantity B 
in (10) remains a fluctuating variable. To eliminate the 
stretching part from the conditionally averaged coeffi- 
cient (A) n — this being critical for explaining the observed 
similarity in the scaling of structure functions between B 
and 9 — one does not need to satisfy conditions (11) and 
(12) for all realizations of the magnetic field B, but only 
for the subset of realizations that gives the main sta- 
tistical contribution to the structure functions (3). Let 
us name this subset of realizations as I. The structure 
functions (3) depend on the statistical properties of the 
increments with respect to t, namely AB T , belonging to 
the inertial range of scales. One of the consequences of 
intermittency is that the statistical properties of the in- 
crements are essentially different from those of the field B 
itself. Therefore, the subset I need not generally coincide 
with the subset G, say, that gives the main statistical 
contribution to the global average {niUjdvk/ dxi) . This 
means, in particular, that the conditions (11) and (12) 
can be valid for the inertial interval (i.e. for subset 7), 
while globally (i.e. for subset G) these conditions could 
well be violated. 



We now use conditions (11) and (12) in the presence of 
the incompressibility condition dvi/dxi = and obtain 

2 >~ < i3 > 

That is, the difference between the passive scalar equa- 
tion (2) and the conditionally averaged equation (10) for 
B is reduced to pure "friction" with the friction coeffi- 
cient given by (13). Equation (10) can then be reduced 
in Lagrangian variables to 

if = w» B > ( 14 ) 

with the "multiplicative noise" (A) n given by Eq. (13). 
Weak diffusion of Lagrangian "particles" can be de- 
scribed as their wandering around the deterministic tra- 
jectories. Introduction of a weak diffusion is equivalent to 
introduction of additional averaging in Eq. (14) over ran- 
dom trajectories [29(. The small parameter r\ in (13) and 
(14) will then determine a slow time in comparison with 
the time scales in the inertial interval and will therefore 
not affect scaling properties of B in the inertial interval. 
This explains the similarity of scaling between B and 6. 

THE MAGNETIC FIELD GRADIENTS 

If one can neglect the last term in the right-had side 
of the equation (10) in the inertial range, then one can 
readily derive equation for the magnitude gradients G = 
VB 

~dt ~ {v ^dx- ~ ~dx~ Gl + V ~dxj> (15) 

The magnitude G of the gradient is determined by G = 
Gg, where g is the unit vector with its direction along 
vector G. Multiplying both sides of Eq. (15) by gi, mak- 
ing summation over i, and taking into account of the fact 
that g? = 1, we obtain 

^ = -((v) n -V)G + 7 ? V 2 G-A'G, (16) 
at 

which is formally similar to Eq. (2) (cf. Eq. (8)) except 
for the last term in (16). The coefficient A' in this term 
has the form (cf. Eq. (9)) 

' / '' / '^- • "UrJ ■ (17) 

One can see remarkable similarity to the situation de- 
scribed in previous section. This similarity suggests ap- 
plying ensemble conditional average similar to that de- 
scribed above. Fix the magnitude G in the vector field 
G = Gn while performing the average over all realiza- 
tions of the direction vector field g permitted by equation 
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(15). Let us denote this ensemble average as (...) g . From 
the definition, this averaging procedure does not affect G 
itself, but modifies the velocity field (v) n , which in turn 
modifies the coefficient A' in Eq. (16). We may write 



% = -(«v> B >. 



V)G + tjV 2 G + (A')gG. 



(18) 



The further analysis can be performed precisely as it 
has been done above for magnitude of magnetic filed it- 
self. Two main consequences of the directional average: 
isotropization and smoothing of the velocity field (v}„ 
and the " nullification" of the production term in the con- 
ditionally averaged equations, can have different propor- 
tions in these two cases. For the above considered case 
with magnitude of magnetic field the conditionally av- 
eraged (on the directions of the magnetic field) velocity 
was still enough strongly fluctuating to belong to the in- 
ertial range paradigm. It is possible, however, that the 
conditional average on the directions of the magnitude 
gradients can smooth the already smoothed velocity field 
(v) n to a substantially non-fluctuating state. In the last 
case the essential point is that the twice conditionally 
averaged velocity ((v) n ) g is smoothed substantially in 
comparison with v, while the fluctuation of G itself is 
still rapid in the diffusion-advection equation (18) (be- 
cause it remains in tact under the conditional average, 
by virtue of its definition) . Under these typical circum- 
stances, the natural expectation (see, for instance, [3l| 
and references therein) is that the space autocorrelation 
function can be characterized by a logarithmic behavior 
[33 | given by 



C(r) 



(G(r)G(Q)) 
<G(0)2) 



In 



L 



(19) 



The result owes itself to the pioneering work of Batchelor 
[3^ . 34 1 who applied this general idea to the viscous- 



convection range of passive scalar fluctuations with large 
Pr. While the two contexts are quite different, they are 
the same in the sense that the velocity field is smooth 
whereas the advected quantity is strongly fluctuating. 

For turbulent flows the Taylor hypothesis is generally 
used to interpret the data. This hypothesis states that 
the intrinsic time dependence of the magnetic field can 
be ignored when the turbulence is convected past the 
probes at nearly constant speed. With this hypothesis, 
the temporal dynamics should reflect the spatial one 
Equation (19) can be rewritten in the temporal terms 
(r -> r L — > t ) 



(G(r)G(O)) /ton 
C[T >- (G(0)2) ln U 



(20) 



This is seen from Figs. 6 and 7 (in the semi-log scales) 
to apply quite precisely for the solar wind data (the solid 
straight line indicates the logarithmic dependence (20)) 
obtained from Advanced Composition Explorer (ACE) 
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FIG. 6: Autocorrelation function C(r) of the magnitude of 
the gradient of B plotted against logr, for the ACE- 1998 
solar wind data. 



O 




FIG. 7: The same as in Fig. 6 but for 2004 year. 



and 



satellite magnetometers for the year 1998 (Fig. 
for the year 2004 (Fig. 7). 

The so-called pumping scale [32| L (or r ) in (19), (20) 
can be defined from the Figs. 6 and 7 as intersection of 
the solid straight line with the temporal r co-ordinate 
axis: tq ~ 24h ± 2h. 



GEOPHYSICAL CONSEQUENCES 

In order to get away from the effects of the Earth's 
magnetic field, the ACE spacecraft orbits at the LI li- 
bration point which is a point of Earth-Sun gravitational 
equilibrium about 1.5 million km from Earth and 148.5 
million km from the Sun. Let us recall that the Earth's 
magnetosphere extends into the vacuum of space from 
approximately 80 to 60,000 kilometers on the side to- 
ward the Sun. Therefore, the observed in Figs. 6,7 ap- 
proximately 24h pumping time-scale of the magnetic field 
gradients cannot be a consequence of the 24h-period of 
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the Earth's rotation. The pumping scale certainly cor- 
responds to a process in the solar wind plasma itself. 
It could be a characteristic scale of the phase coherent 
structures in the solar wind, which can be distinguished 
from incoherent fluctuations in this case (see more 
about characteristic periods in solar wind plasma in the 
recent review [13] )• It is now believed that energy of these 
coherent (pumping) large-scale structures is released via 
shear instabilities in the solar wind plasma. 

It is possible that for the very long time of the 
coexistence of Earth and of the solar wind the weak 
interaction between the solar wind and Earth could 
lead to stochastic synchronization between the Earth's 
rotation and the characteristic time-scale (the pumping 
scale) of the solar wind magnetic field gradients. This 
synchronization could transform an original period of 
the Earth's rotation to the period close to the pumping 
scale, which we observe at present time. 
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T. Gotoh for providing Ref. [lij before publication. I also 
thank K.R. Sreenivasan for inspiring cooperation, and D. 
Donzis, J. Schumacher, and V. Steinberg for comments 
and suggestions. 
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